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1. Introduction 


The use of formal series as a tool for stud ying inte grators of ordinary differ¬ 
ential equations has a long history, see e.g. SSMlhl ] for an overview. Butcher 
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But72l | introduced series indexed by rooted trees, B-series, to study Runge- 
Kutta schemes and generalization s of these. The algebraic struc ture of B - 
series has been thoroughly studied |CHVlCll : lMur99l | and others. In |MKV15I |. 
a generalization of B-series, called aromatic B-series, was introduced as a 
classification of numerical methods that are affine equivariant, that is not 
affected by linear changes of coordinates. This paper derives counterparts 
to some of the results concerning algebraic structure of B-series to aromatic 
B-series. 

In particular: 

• The composition law of B-series CHVlfll . Theorem 3.1] is mirrored by 
the composition of aromatic S-series theorem 14.11 

• The relationship between S-series and B-series Mur99l . Theorem 3] is 
mirrored by the slightly more complicated relation in theorem 14.21 

• The substitution law for B-series [CHVld . Theorem 3.2] is mirrored by 
the substitution law for aromatic B/S-series theorem 15.31 


2. Aromatic B-series 


In MKV15| . aromatic B-series are introduced as a tool to describe numerical 
integrators which are equivariant under invertible affine maps M” —M”. In 
this section, we define aromatic B-series and aromatic S-series. 

Aromatic B-series are series indexed by equivalence classes of certain di¬ 
rected graphs. For our purposes, a directed graph 7 = (R, E) is dehned by 
a finite set of vertices or nodes V, and a set of edges E C V x V . We say 
that the edge {vi,V 2 ) goes out of vi and into V 2 - A subgraph of 7 , is another 
directed graph (IT, F) where IT C R, F C IT x IF fl i?. In this definition of 
graph, we allow the empty graph with 0 vertices, and self-loops. 

On a given graph, we define the direct predecessor function from R to 
the power set of R by ui G TT-y{v 2 ) iff {vi,V 2 ) G E. 

Two graphs are equivalent and we write 71 = (Ri,Fi) ~ (R 2 ,T^ 2 ) = 72 if 
there exists a bijection g: Vi —>■ R 2 such that {g x g){Ei) = F 2 . 

The automorphism group of a directed graph is the set of permutations 
( 7 : R — 7 > R satisfying that {g x g){E) = E. 

Our object of study is the equivalence classes of directed graphs satisfying 
an additional criterion. 


Definition 2.1. An aromatic forest is an equivalence class of directed graphs 
where each node has at most one outgoing edge. We denote the set of aro¬ 
matic forests as AE. A root of an aromatic forest is a node with zero outgoing 
edges. The set of roots of the aromatic forest cj) is denoted r{(f)). 

We will occasionally make statements where an aromatic forest is a graph 
(as opposed to an equivalence class of graphs). Take these statements to mean 
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a member of the equivalence class, the independence of choice of member is 
either obvious or explicitly stated. 

It follows from the definition that an aromatic forest consists of connected 
components, each of which has either (a) one root, in which case the con¬ 
nected component is a rooted tree, or (b) no roots, in which case it contains 
exactly one cycle. 

Definition 2.2. We define the following subsets of AJ-'. 

Ji,... }, the set of aromatic forests with no roots. 

• A' = Ji,... }, the set of connected aromatic forests with no 

roots. 

• AT = •, I, (5)», V..., the set of aromatic forests with exactly one root. 

• T = •,t,V..., the set of rooted trees, or connected aromatic forests 
with exactly one root. 

• T, the set of loopless aromatic forests, which are (unordered) multisets 
of rooted trees. The set T is also called forests. 

• Tk, the set of forests with exactly k roots. 

We note that AT = Ax T. 

Definition 2.3. An aromatic forest f represents a elemental differential op¬ 
erator, Ff{(j)) acting on smooth functions defined in the following manner: Let 
(j) = {V, E) For each node p & V, form the factor ^, where T(p) = iqiiq 2 ''' 
is the multiindex defined by the direct predecessor function on f, that is 
qi,q 2 , ■ ■ ■ are the vertices such {qj,p) G E The upper index on / corresponds 
to the vector components of / and the lower are partial derivatives with re¬ 
spect to the coordinate directions, j■ ■ ■ dxi^^. Then, 

for each root s, form the factor di^ = d/dxi^. Finally, multiply the factors 
and sum over repeated indices. 

We can write 

w = n f'T, n 9- 

p&V{y) s&{'y) 

where we sum over repeated indices. 

Example 2.4 (Elemental differential operator). Let 7 be the tree with 
indices 

i 3 k 

I 

n 

Then FfiP = f Pf ffd^. 
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For (p a tree, the elemental differental operator Fj[(j)) corresponds to the 
elemental vector field from the B-series theory. 

In general, an elemental differential operator Ff{(j)) is a differential op¬ 
erator of the same degree as the number of roots of p. Thus, if 7 € 
Ff{'y) is a scalar field and p = 7 rir 2 ■ ■ - Tm, where ■'y ^ A, Ti & T, then 
FfiP) = Ff{'y)Ff(Ti) ■ ■ ■ Ffijm) is the product of the scalar field Ff{'y) with 
the (frozen) composition of m vector fields Fj(ri),... Ff{Tm)- 


Definition 2.5. In MKVlf)l | an aromatic B-series is defined as a series 


indexed by aromatic trees. It represents the vector field 


Bf(b) 


E'-"' 

t€AT 


(T(r) 


Fy{r). 


The normalization constant a is here defined, for an aromatic forest (/>, to be 
the size of G^, the graph automorphism group of p. |r| denotes the size or 
number of vertices in r. 


An aromatic B-series is a series expansion of the modified vector field of a 
numerical integrator. As for B-series, it is possible to show that the numerical 
integrator can also be expanded in an aromatic B-series, y y + Bf{a)[y]. 

Proposition 2.6. The differential equation 

y ^ E Pp ri{T){y), y(0} 

rtir 


has solution 


y{t) = yo+ 


for some o: k x AT k, 


r,=AT 

a: (t,r) 1 -^ at(r). 


Proof. Let R{t) = y{t) — yo- By integration and Taylor’s theorem, we see 
that 


^ ' fc =0 


r&AT 


Ljo 


R{s)^^ds 


where Ff{T)^^\yo) is the A:’th derivative of at yo taken as a function 

^E. Since Ff{TY^1{Ff{9i)^ ...0Ff{dk)) is again a sum of elemental 
differential operators of the form Ff{Q), where |0| > |r|, we can write R{t) 


as a 


sum -recursion on |r|. 


□ 


To be able to compose aromatic B-series, we want to describe the effect of 
such an integrator to a function, i.e. evaluate G(y -|- Bf{a){y)) for smooth 
G. To do this, we introduce aromatic S-series, which mirrors the S-series of 
Murua Mur99 |. 
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Definition 2.7. Let o: AT —?> k be a series indexed by aromatic forest. Let 
the aromatic S-series of a be the operator acting on functions defined by 

s/WPi = E 

Calculating the aromatic S-series corresponding to the mapping y ^ y + 
Bf{a){y) requires some algebraic machinery which will be developed in the 
next section. 

For now, we define the composition product of two S-series. 

Definition 2.8. For a, b: AT k, let their composition product o • 6 be the 
aromatic B-series defined by Sf{a ■ b)[G] = S'/(6)[5/(a)[G]]. 

An important building block of the composition product is the composition 
of two elemental differential operators. 

Definition 2.9. The composition of two aromatic forests 4>iO(p2 is the formal 
sum of aromatic forests defined as follows: Let * be a dummy node not 
contained in either or V{<f) 2 ) and {V{<f) 2 ) the set of functions 

p: ^ V{4>2) U -k. Let 

4>1°4>2= 

peiVi<i>2)u*Y(M 

where ‘h((/)i, (j) 2 , p) is the aromatic forest with nodes V ((/>i) U V{ 412 ) and edges 
F;((/>i) U E{42) U Ep, where Ep = {{p, p{p)) s.t. p € r(^i), /9(p) / *}. 

Informally, the composition product is formed by the sum over all possible 
ways to add edges from some (possibly none) roots of 41 to the nodes of 42 - 

Lemma 2.10. 

Ef{4i)[Ef{42)[G]]=Ef{4io42)[G] 

where the right hand should be interpreted as if Ff is a linear function. 
Proof. By 12.31 T/Ida 1 is a sum of products of the form 

n n 

peV{(l>i) ser(0i) 

and similar for 42- Use the product rule to expand 

i^/(*)[f/(*i)[Gii = n s'T n n n o.,. 

The result is a sum of terms, where in each term, each of the operators di^ 
is applied to (a) one of the factors ^ or (b) one of the factors Gi^. 

Referring to the definition of the elemental differential operators and the 
composition, case (a) corresponds to adding an edge s ^ q, while (b) corre¬ 
sponds to not adding an edge from s. 

The sum is then taken over all possible choices. □ 
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3. Coalgebra structure 


We will be interested in series indexed by aromatic forests, and composition 
of such series. A convenient environment for studying series indexed by 
elements of a certain set S is the dual space of the free vector space over the 
set S. 

Let 


^ ^ 7 ) Hjr = © kw and C XF = ® Hjr = kt/i, 

't&A 4>eAF 

denote the free vector spaces over the sets A, T and AT. Hjr is as a vector 
space isomorphic to the Hopf algebra of rooted trees, and (7^ is the vector 
space of commutative polynomials in the variables A!. 

In the following, we will describe coalgebrE0 structures on Cx Hjr and 
CxF- We recall that a coalgebra is a vector space C equipped wi th a co - 
product A: (7 —)• C ( 8 i <7 and a co-unit e: (7 —>■ k. Following Sweedler (Swe69|, 
the coproduct is written A(c) = i 8 > C( 2 ). For a precise definition of a 

coalgebra, see Appendix iBl or |Swe 6 fl |. 

The coproduct in (7^ is a scaled coproduct of multiplication of commu¬ 
tative series in A'. Let a € be identified with the formal infinite series 
'^■yeA 71^7 correspondingly for b E C^. Let a ■ b G be defined by 


o • 6 ( 7 ) 

2^^77W7 = 


76.4 


a{j) 


E?M71 


71 


(t(7i 




72 


o-(72 


where the right hand side is the formal multiplication of series. 

It is clear that a ■ 6 ( 7 ) is a hnite sum of terms of the form a( 7 (i))K 7 ( 2 ))) 
and we define the coproduct Ayi( 7 ) = J2{-y) 7 (i) ® 7 ( 2 ) by 

a-b{-r) = ^a(7(i))6(7(2)). 

(7) 

Co-associativity of A _4 follows from associativity of formal commutative mul¬ 
tiplication of series. Also, we define the counit e _4 to be the map 1*: ( 7_4 —)■ k. 

The copro duct of Hjr is the dualization of the Grossman-Larson product 
as dehned in CK98]. 

Both C_A and Hjr are graded connected coalgebras where the grading is 
induced by the number of vertices of the graphs, and {C_a)o = [Hjr)^ = kl ~ 

k. 

Since Cj,jr = C_a G Hjr, we can identify 


^All (co)algebras considered in this paper are (co)associative and (co)unital, and we will omit the 
qualihers. 
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where CX) denotes linear functions from Hjr to in the usual 

way: a: Hjr —)• k is mapped to a'\ Hjr by 

{a'{4>),c) = a{c® (j)). ( 1 ) 

Hjr can be equipped with the concatenation product, while is the dual 
of a coalgebra, which is an algebra in a natural way: The product on C'^ is 
defined as 

(o- 6 ,c) = (o( 8 ) 6 , A^(c)) = ^(a,C(i))( 6 ,C( 2 )). ( 2 ) 

(c) 

The unit of is given by the co-unit e^. 

In other words, series indexed by aromatic forests can be identified with 
linear maps Hjr where both sides are algebras. We will see later that 

series which correspond to algebra morphisms Hjr play a special role. 

Definition 3.1. An admissible partition of an aromatic forest 0 is a 
partition of the graph (p into two (possible empty) subgraphs R{p,f)) and 
P*{p<f>)j such that no edges in (p goes from Ript/,) to P*{p^). 

The graded coproduct is defined on aromatic forests as 

</>(!) «) p(2) = P*{p^) ® R{p<p) (3) 

(4.) 

and extended linearly to 

The co-unit is defined by e{p) = 1 if (/> = 1, 0 otherwise. 

Theorem 3.2. (C^, A_ 4 j-, e) is a graded eo-algebra, 

Proof. That A_ 4 _ 7 r is co-associative follows from the observation that (A ® 
id) o A and (id ( 8 ) A) o A both can be interpreted as the splitting of p in three 
parts pi,p 2 , p 3 such that no edges goes from to p 2 nor from p 2 to pi. The 
co-unital property follows from 

e{P*{P4))PiP4) = 

p<i> 

and 

Y P*(P^HPiP^)) = = (t>- 

p<i> 

Finally, the grading is induced by the number of nodes of p. □ 

We note that CXF = Ca ® Rf as a vector space, but A_ 4 j- is not the 
coproduct inherited from this structure, (id 0 r (8 id) o (A^ 0 Aj-), where 
r: C_A 0 Hjr -A Hjr 0 denotes the twist map c® h ^ h ® c. 

We can even equip CxF with the concatenation product U to obtain a Hopf 
algebra. However, as will be shown in the end of section 01 this structure is 
not as closely connected to aromatic B-series integrators as the Hopf algebra 
structure is to B-series integrators. 
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Theorem 3.3. U, 1, A_ 4 jr, e), is a graded, connected bialgebra, thus 

also a Hopf algebra. 

Proof. The compatibility conditions are straight-forward to check, and the 
existence of an antipode follows from theorem lB.il □ 

4. The composition law 

We are now ready to tackle the composition of aromatic B-series. 

Theorem 4.1. Let a,b: ^ h, then Sf{b)[Sf{a)[G]] = Sf{b-a)[G], where 

b ■ a{f) = (6 ® a)(A^(c)) = ^ b{P*{p^))a{R{p^)). (4) 

W 


Proof. By definition 

S,(6)|S;{a)[G]] = ^ 


4>2^-AJ^ 

Using Lemma l2.ini we rewrite 




4>2^AJ^ 




S/((.)[S/(o)[0]] = 


(j)2^-AT 


a{(j)2y 


Y. ^F/(*)|G 1 










a{(j)2)(y{yiy 


EE E 

4>2&AT <i>igAT p£(v(4>i)u*Y‘-‘I‘2) 


b{(t)2)a{4>i) 

a{fi)a{(t)2) 


Ffm(l22Ai,p))[G] 


- E E E 

<i>ieAT<j>2eATY&AT 

where Mi{(fi,(j) 2 ,y) is the multiplicity of fj in 

4>2o4>i= ^ <^{cj)2,4>i,p)- 

peiVi4>i)u*)A<i>2) 

The claim is that [5] is equal to 

b ■ a{y>) 


Sf{b-a)[G]= ^ 


p&AT 


a{y) 


-FjmG] 




'(p^AT (i/)) 


a{y) 


E E E ^-Fyy^K4,,)amF,ma\, 


(f)2^AJ^ '(p^AT 


a{y) 


(5) 


(6) 




















where M 2 {(f)i, < 1 ) 2 , is the multiplicity of (pi 0 (p 2 in Aciip)- 

The theorem holds if 

Mi{(pi,(p2,^) _ M2{(pl,(p2,1p) 
a{(pi)a{(p 2 ) cr(V') 

for all (pi,(p2T'>P- We proceed by proving this equality. 

(1) Ml = 0 M 2 = 0; If edges can be added from (p 2 to (pi to form ip, 
the splitting of ip into (pi (p 2 is an admissible partition and vice versa. 

(2) If Mi,M 2 / 0, cardinality 

of subgroups of, respectively, x and G^ and we will prove that these 
subgroups are isomorphic. 

For the remainder of the proof, let (pi,(p 2 ,'ip, p',p'^ be such that 
^{(p 2 ,(pi,p') ='ip and P*{p'^) ® Rpp'^) = (pi® (p 2 . 

Let G^^ X G (^2 act on {V{(pi) from the left by setting 

{ 91 , 92 )- p{p) = 91 °P° 92^ (p) 
where we define gi{*) = *. 

Now, ^{(p2,(pi,p') - ^{(p2,4’i,{9i,92) ■ p') as graphs, and Mi{(pi,(p 2 ,ip) is 
the size of the orbit of p' under the action of G^^ x The orbit-isotropy 
theorem yields that isotropy group of p' under 

the same action. We call this group Hi. 

The admissible partition of V’j Pi/” fally dehned by the set of vertices 
forming the subgraph R{p^). We let G^ act on the set of admissible par¬ 
titions of Ip by its action on the subgraph R{ip). (That is, by considering 
the simultaneous action on the set of vertices R{ip).) Since G^ acts on 
Ip by graph isomorphisms, the restriction of the action to a subgraph is a 
graph isomorphism of the subgraph. Therefore, R{g ■ p'^) ~ R{p'ip) and 
P*{g ■ p'.^) ^ P*{p'ip) as graphs, and M 2 {(pi,(p 2 ,ip) is the size of the orbit of 

under the action of G. The orbit-isotropy theorem yields that M 2 ( 0 it ^2 p) 
is the size of the isotropy subgroup of under the same action. We call this 
group H 2 . 

We now prove that Hi ~ H 2 as groups. Let ( 51 , 32 ) G Hi i.e. gi o p' = 
p' ° 92 \r{(i> 2 )- Now, define the action of ( 5 'i, 5 ' 2 ) on ip = ^{(p 2 ,(pi, p') by letting 
5i, <72 act on the subgraphs <piT(p 2 - This is a graph automorphism of ip because 
for (r,/j'(r)) € {g 2 {r),gi o p'{r)) = {g 2 {r), p' o g 2 {r)) € E'^. It is in H 2 

because ( 51 , 32 ) niaps i?(p^) (pi to i?(p(^). 

The converse isomorphism H 2 —> Hi follows from considering g G H 2 and 
doing the appropriate restrictions to get ( 51 , 52 ) £ G(j(,j x G^j- For {p,p{p)) G 
E{ip) , we get (52(p),5i o p{p)) G E{ip). As g 2 {p) G (p 2 , 5i o p{p) G (pi, we 
must have 51 o p{p) = pog2{p). Therefore ( 51 , 52 ) G Hi. 
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As Hi ~ i/ 2 ) we get 


a{^i)a{^2) ^ 1^ I 1^ I 

Mi{(j)i,(l>2,'lp) ^ ^ M2((/>i,02,V’)’ 


□ 


Theorem O says that the composition product (8i C^jr C'xF 
the dual of the coproduct in C^F- As the dual space of a graded connected 
coalgebra, the invertible elements of C’^ are {a € s.t. a(l) / 0}, and 
the inverse is given recursively by 


a (1) = l/a(l), a (</>) = - 


«(</>) +Ew «(</>')« i(p") 

a{l) 


where ® 4>" = ^{'P) — — <^<8)1 is the reduced coproduct. 

We are now ready to calculate the aromatic S-series oi y i-» y + Bf{a)\y]. 
The proof closely follows the similar proof for S-series in Mur99j | . 


Theorem 4.2. Extend a: Cj\j- —> k to a function a: Cjijr h by requiring 
that a and a coincide on AT, and that a': Cjr —>■ defined as in ([T]), is an 
algebra morphism. Then G{y + Bf{a)\y\) = Sf{d)[G]{y). 

Proof. AT = A X T, so we can rewrite the sum 


B,[a) 


E 

t£AT 


a{T) 

a{T) 


= E E 

esT 7eT 


a{l0) 

cj( 76 ») 


Ff{l0), 


Now, we make two observations: (i) Fj{'y9) = Ff{^)Ff{6). (ii) The auto¬ 
morphism group of 70 cannot mix connected components with 0 and 1 root, 
so cr{'y6) = a{'y)a{6). Collecting terms, we can write 


Ffia) = 

9eT 


E 

76.4 


C7(7) 


T/(7) 


1 

a{e) 


Ff{ 0 ). 


We now write a{'y6) = {a'{0),^) where a'{0) G is as in ([T]). With a slight 
abuse of notation, we write. 


Sf(a'm 


E 

7eT 


{a'{0)T) 

cj(7) 


Ffil)- 


Sf{a'{9)) is a scalar function of y. We can now write 


B/(»)= E 

e&T 


Sfia'ie)) 

a(9) 


Ff{9). 
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We see that this is simply the expression for a regular B-series, with the 
constant coefficient a{9) replaced by a y-dependent scalar value Sf{a'{9)). 
Using Taylor’s theorem, we get 


G{y + Bf{a){y)) =G [y + Bf{a){y)] . 




k=l 


The multilinearity of G^^\y) allows us to rewrite the above equation. We 
suppress the dependence on y 


G{y + Bf{a){y)) - G{y) = 




In the above equation, ui • • • is a fc-tuple of rooted trees, u is the number 
of distinct trees among ui,...,Uk, and yi,..., are the multiplicities of 
each distinct tree. 

If (/> = tti • • • Um, Gfj, is the semidirect product of permuting identical trees 
among ui,...,Um and Gui x ... x Gu^ , so 


a{(j)) = yi! • • • yjaiui) ■ ■ ■ a{uk). 
We define a': Hjr —>■ to be an algebra morphism by 

= a'{ui) ■ ■ ■ a'{uk) 


where the product is taken in G^. For 4> = 1 , this simplifies to the empty 
product a'(l) = e^. 

For any 71,72 G A, we have ^^( 7172 ) = Ff{'yi)Ff{'y 2 ). Consequently, we 
have 5j-(a'((/>)) = Sf{a'{ui)) ■ ■ ■ Sf{a'{uk)), and 



□ 


As for B-series, we can recover Bf{a) = Sf{a)[I] — L 
The aromatic S-series a appearing in theorem 14.21 are, in general, not char¬ 
acters of (Cyij-, U) with values in k. To see this, consider the aromatic graph 
The requirement that a': Gjr —>■ is an algebra morphism implies 


o(s5) = (a'(l), (5) = (e^, = 0, 
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on the other hand, a((5)») = a((5)») can be nonzero. 

The S-series in theorem 14.21 are characters of Hjr with values in C^, and 
such characters form a group when equipped with the convolution product 

a' • b\^ 2 ))- 

(t) 

However, the composition product of aromatic S-series theorem 14.11 is a dif¬ 
ferent product. 


5. Substitution law 


In GHVlOj, the effect of replacing the vector field / in an integrator with 
a different vector field /, expressed as a B-series, is studied. The same 
procedure can be applied to aromatic B-series. In this section we express the 
substitution law for aromatic B-series. 


Definition 5.1. A partition p‘^ of an aromatic forest (f> is a partition of the 
graph (j) into subgraphs P{p^) = {9\, ..., Om} such that 0i,..., 0m € ■AT. 
We view P(p'^) as a multiset of aromatic trees. The set of partitions of </> is 
denoted V{(j)) The skeleton xip'^) of ^ partition is the aromatic forest which 
is obtained by collapsing each of the subgraphs in P(p'^) to a single node. 

Constructively, a partition is defined by cutting a subset of the edges in (/>, 
obtaining a set of connected components, and then assigning each rootless 
connected component, to a connected component with root (i.e. a rooted 
tree). 

Definition 5.2. For the purposes of the proof, we also define a labelled 
partition pT as a partition where we assign each subgraph 0 i,... to a vertex 
of V{x{p*(t>)). For a labelled partition, we let P{p*^) = {Oi,... ,0m) be a 
sequence of aromatic trees. The set of labelled partitions of (p is written as 
P*{cP). 

For a partition p^, the number of possible labellings is cr(x(p'^)). 

Theorem 5.3. Let b: AT h be an aromatic B-series and a: AT h be 
an aromatie S-series, then SBf(b){cL) can be expressed as an aromatic B-series 
in f, SBf(b){cL) = Sf{b-ka) with coeffieients 


b-ka{(p)= «(x(/)) n 


Proof. We prove the theorem for 0 = 7 * where 7 is an arbitrary graph. The 
full theorem follows from linearity in a. 
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By the definition of the elementary vector fields, Sg{a) = where 

Pg{l) = nfcGr(7) ^k, 

Now, 


9eAT 


so 


r,w=EE-E n n 


cr(9j) 

6 l O 2 dm jGy{"j) ^ 


-^U) 


k€r{'y) 


In the above equation, njGy( 7 )^ j is a linear combination of differ¬ 
ential operators. The aromatic graphs corresponding to these operators are 
formed from the graphs 61 ,... , 6 m in the following manner: For each edge 
ei = {u,v) in E{'y), choose a vertex in 6 ^ and add an edge from the root of 
9u to that vertex. The sum is then taken over all possible choices of vertices. 

We can denote the aromatic forest formed by a particular choice of vertices 
by 

r(6»i, ...,9m-,ji,-- •, jn) = r(0; j). 


where ■ ■ ■ jn are the vertices chosen for edges ei,... ,en S Recall 

that for ei = {u,v), we are restricted to choosing ji in 9^. 

With this notation. 


Sg{a) 


1 

cj(7) 


E E 

eeAT^ j 


Ui K^i) 


^/(r(0;j)), 


where the sum goes over all possible choices of vertices. 

We now collect equal terms. Collecting j for which r(0; j) = (j) are equal, 
we get 

( 8 ) 

Where Mi ( 7 , 0, (p) counts the number of j such that r(0; j) = cp. 

What we seek to prove is that Sg{a) is equal to 

= E n wow. 

0 ^ P^eVi(j>) ees{p^) 

X(p‘^)=7 


We rewrite this by extending the sum over all labellings of the partition, 
and dividing by £ 1 ( 7 ) to compensate. 


Sfib^a) 


E 


1 

(T{(p)a{j) 


p*<p£P*{(f)) eeS{p*'^) 
x(p*‘^)=7 
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Collecting terms for which P{p*'^) = (0i,..., 9m) = ©, we get 


where M 2 {'y,Q, 4 >) counts the number of labelled partitions p*^ such that 
S{p*<t>) = 0. 

Comparing ([5D with Q, Sg{a) = Sf{b-ka) provided that 
M2(7,0,(/)) ^ Mi(7,0» 

We proceed in a similar manner to the proof of theorem 14.11 

If there is no way of adding edges to the graphs 6 i,... ,6m that results in 
(j), there is also no way to obtain 9i,... ,9m by cutting edges in cj) and vice 
versa. In this case, both sides of the above equation is zero. 

If non-zero, we can identify each side of this equation with the reciprocal 
of the cardinality of an isotropy subgroup. 

(i) Gfj, acts on labeled partitions of </> in a natural manner. Let p*^ be such 
that xiP*^) = 7 ) and S{p*^) = {9i,..., 9m)- The action of conserves both 
x{p*^) and S{p*^) = {9i,... ,9m), and all partitions satisfying S{p*(j)) = 0 
is obtained in this manner. Therefore, M 2 ( 7 , ©,</>) is the size of the orbit of 
p*’f’ under the G^ action. 

By the orbit-isotropy theorem, is the size of the isotropy subgroup 

of p*'^, that is: the subgroup of G^ which does not mix the vertex subsets 
corresponding to 0i,..., 9m, and in addition, within each vertex subset V(9i), 
any vertex whose in-degree is higher in (p than in 9i, has to be fixed under 
the subgroup. 

(ii) j denotes a sequence of vertices in 9i,...,9m with some restrictions, 
X... X Gg^ acts on the set of such sequences by moving the vertices within 

each subgraph. Let j be such that r(0; j) = 4> The action of Gq.^ x ... x 
preserves r( 0 ; j), and all sequences satisfying the requirement are obtained 
in this manner. Therefore Mi( 7 , 0,0) is the size of the orbit of j under the 
action of Gq^ x ... x Gq^ . 

By the orbit-isotropy theorem, is the size of the isotropy group 

of j under this action, that is the subgroup of Gg.^ x ... x Gg^ which, in each 
component preserves j n V{9i). 

Since the vertices in j are exactly those which have higher in-degree in 0 
than in 6i, the two subgroups described in (i) and (ii) are isomorphic. 

□ 

Proposition 5.4. (a) (G^-y-,*) forms a monoid with identity** which acts 

on G^jr by algebra morphisms, in particular, the following identities 
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hold for a, ai, 02 € C*^jr, b,c & ^AT’ <^15 € k. 

b -k [aiai + aiai) = aibk ai + a 2 & * 
ck {bk a) = (ckb) k a 
bk,* = b 

★ (X = tt 

bk (ai ■ 02) = {bk ai) ■ {b k 02) 
bka~^ = {bka)~^ 


( 10 ) 


(b) The set {b G C'^t- s.t. b{») / 0} with operation k forms a group. 

Proof (a) All of the formulas follows from Sf{bka) = SBf(b){o-)i and Sf{a- 
h) = Sf{b)oSf{a). 

(b) From (a) we know that •* is both a right and left identity on The 

right inverse of an element b G b(») ^ 0 is given by 


= Hx(p*)) n ''"‘W 

e&s{p'i>) 

^trivial 

where in the last sum, the trivial partition with x(p'^) = •, S{p'^) = r is 
omitted. All elements in {b G Cf^j- s.t. 6 (») 7 ^ 0} allow a right inverse 
in the same set, therefore a right inverse is also a left inverse. 

□ 


Divergence of aromatic B-series 


A special case of the substitution law is the divergence of an aromatic B- 
series. This is obtained by setting a = ( 5 * in theorem 15.31 and we set Vb = 
bk {!)*. Since Ff{f>) = div(/), Sf{b k {!>*) = div Bf{b). 

In this case x(p'^) is required to be the 1-loop (5, and P{p'^) can only consist 
of a single aromatic tree. V 6 is then only non-zero on aromatic graphs with 
zero roots, and we get 

Vb{cf)= Y. b{f,), 

v£V{(l>) 


where cfv is the aromatic tree obtained by deleting the e dge going out of v. 
For B-series, this formula is considered in CMOTI : llQT07 |. 

Now, assume that / is a divergence-free vector field div / = Ff{it>) = 0, and 
that we apply an aromatic B-series method (i.e. an affine-equivariant, local 
numerical method) to the differential equation y = f{y). The modified vector 
field given by Bf{b), and the method is volume-preserving if V 6 = bk ^ 
non-zero only on </> G AF that have (5 as a subgraph. 


IS 
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Affine-equivariant pseudo-volumepreserving 
integrators 


Volume preserving methods are well known, see for instance FS95I : lQui95l |. 
but all such methods rely on some partition of the vector field along coordi¬ 
nate axes, and are therefore not affine equivariant, and not aromatic B-series 
methods. 

It is known CMOTI : llOTO?! that B-series methods (apart from the exact 
integra tor) cann ot be volume-preserving, but that aromatic B-series methods 
can be MKV15l |. However, it has so far not been possible to find an equation 
defining the update map for a volume-preserving aromatic B-series method. 

With the computational tools developed in this paper, one can attempt 
to find methods that are volume-preserving up to a certain order by using a 
pre-processing vector fields. 

An example of a method which is second order, but preserves volume to 
fourth order is detailed below. 


Example 5.5. We base our method on the implicit midpoint method, whose 
modified vector field is Bf{a), where 


^ I 

CL — • -\- 


1 

12 


r 


— V* I -|- remainder. 


Here the remainder contains terms of order 5 and higher. 
And define a pre-processed vector field 


F{y) = f{y) + 


(1 


12 V 2 


■tr(/'( 2 /)^)-/'(y)^ /(y). 


In the notation of aromatic B-series, F = Bf{b) where 


b 




The integrator is defined by substituting F for / in the implicit midpoint 
method 

I u T? ( F Vn+l ^ 

yn+i = yn + hF[ - - - 1 . ( 11 ) 

Implementing this method requires evaluating the Jacobian f'{y) and solving 
a non-linear equation. 

The modified vector field of m is Bf{b-ka) and the divergence of the 
modified vector field is Bf{b-ka-kil>*), which we can calculate for the smallest 
rootless graphs. Using © 


b-k a 




V*) -|- remainder. 
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which we can reuse in (6 * a) * (5* 

1 

b-k a-k ^* = C>* — —P + remainder 

where the remainders contain terms of 5th order or higher. 

If / is divergence free, Ff sends any aromatic forest containing the sub¬ 
graph i?) to 0 , so we see that while Bf{bk a) agrees to / to second order, 
div Bf{bka) = Sf{bka-k^*) disappears to fourth order. The integrator (fTT]l . 
while second order, conserves volume to fourth order. 

The integrator in example 1 5.5 1 is. to the author’s knowledge, the first exam¬ 
ple of a pseudo-volumepreserving affine-equivariant integrator presented in 
literature. The integrator was derived using the substitution law for aromatic 
vector fields theorem 15.31 Further calculations using the substitution law, not 
included in the paper, have unfortunately shown that arriving at a true vol¬ 
umepreserving affine-equivariant integrator is unlikely to be achievable by a 
similar process. 
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A. Tables for composition and substitution laws 


(p 

^at{4>) 

1 

1 (g) 1 

• 

1 0 • “h • 0 1 


l0S5 + i50l 

: 

l0j“h*0*H“J0l 




1 (g) (^5 + (g) 1 

V 

1 0V “1“ 2»0j-t-#»0»-|-V 0 1 

i 

1 0iH“*0lH“l0*H“i0 1 


l(g)Y + 2«(g)if) + »«(g(!?)+Y<gl 

1 

l(g)|+.(g)(?)+l(g!5+i(g)l 

h 

l(gB + »<8)i^^ + B<8)l 


1 (g (g 1 

•V 

1 0 -1-3«0V -|-3»«0l-|-»»»0»-|- 0 1 

V 

• 0V H-»0iH-l0lH-»»0lH-I»0»H-^ 0 1 

Y 

i(gY + 2 .gi+,.gi+vg.+Ygi 

• • 

i 

l0lH“»0iH“l0lH“i0»H“l0l 

• 

1 0 “1“3#0^ -l-3##0i?) -l-##*0^ -\- 0 1 

• • 

k 

• 0^ + *0^+l0^+**0^+t*0^>+K0l 

% 

l(gX + 2.g|)+..gif) + Vg(5+Jgl 

I j 

1 

lgi^+.g6+tg(?)+ig!5 + (5gl 


i(gB + 2.gB + **gi^ + Bgi 

• • 


igB+«gB+ig(» + Bgi 

B 

lgS + 2.g!b + **gi^ + Sgl 


igB*> + «g(***) + B*>gi 




Table 1: The coproduct ^at 
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: 

(5. 

V 

i 

S5l 

h. 

(5(5. 

(5V 

(5i 

(V)I 


6 * a(r) 

«(•)&(•) 

a(.)6(J) + a(.)6(.)' 
a(.)6((5.) + a((5.)6(.)2 
a(.)6(V) + 2a(I)6(.)6(I) + a(V)6(.)=' 
a(.)6(i) + 2a(I)6(I)6(.) + a(i)6(.)=' 
a(*)6((5l) + 2a(I)6((5*)6(») + a(l)6(l)6(*) + a((5l)6(»)^ 
a(*)6({55») + 2a((5»)6(I)6(*) + a((55»)6(»)^ 

a(*)6((?i») + a(I)6((5*)6(») + a((5»)6((5»)6(») 
+a((5.)6(I)6(.) + a{h»)b{»Y 

a(»)6((5(5*) + 4a((5*)6((5*)6(») + a((5(5»)6(*)^ 

a(.)6((5V) + a((5.)6(V)6(.) + 2a(I)6((5.)6(I) + 2a(I)6((5l)6(.) 
+a((5l)6(I)6(.)2 + 3a(V)6((5.)6(.)2 + a((5V)6(.)^ 

a(«)6((5i) + a((5«)6(i)6(«) + 2a(X)6((5«)6(t) + 2a(I)6((5t)6(«) 
+2a((5l)6(I)6(.)2 + 3a(i)6((5.)&(.)' + a((5i)6(.)^ 

a(»)6((55l) + 2a(I)6(^*)6(») + 2a((5*)6(I)^ 
+2a((5l)6(I)6(.)2 + a((V).)6(I)6(.)2 + a{(^l)b{,f 

a(»)b{^i) + a(I)6((5l)6(») + a((5»)6((5»)6(t) + a((5»)6(I)6(») 
+2a(I)6((?..)6(.) + a(V)6((5.)6(.)2 + a(i)6((5.)&(.)' 

+a((5l)6((5»)6(*)^ + a(^*)6(I)6(*)^ + a(^»)6(I)6(»)^ + a{^t)b{») 


4 


Table 2: Substitution law (1 root) 
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A 

T 


S5 

h 

I 

h 

hc> 


b -k a(A) 
a(l) 
a{^)b{,) 

a(s5)6(t) + a(s5)6(i5») + a(s?>)6(*)^ 

2a(s5)6(I) + a(sJ^)6(.)2 
2a({5)6(s5*) + a(s5s5)6(»)^ 

2a{())b{hf>) + a{())h{*f) + 2a{^)b{()»)b{») 
+2a{h)b{l)h{.) + a($)h{.f 

a(s5)6(^*) + a(s5)6({5t) + a(s5)6(i) + a{^)b{^»)b{*) 

+a((?))6((5«)6(«) + 2a((?))6(I)6(.) + a(|))6(.)^ 

a(i5)6(i) + a(<5)6(V) + a{C>)b{^.) 
+2a((?))6(.)6(t) + a((J5)6(.)6(I) + a(^)6(.)^ 

3a(!5)6(i) + 3a((J5)6(I)6(.) + a(!i?5)6(.)^ 

a{i!))b{^») + 2a((5)6((5t) + 2a((5(!?))6(t)6(») 
+2a(sJ^)6(s5.)6(.) + a({r5i5)6(.)3 

a(s5)6({5{5») + a({5)6({5l) + a(i5)6(s?>.) + 2a(s?>)6(s5*)6(») 
+a((5(5)6(I)6(.) + a(if))6((5.)6(.) + a((5!5)6((5.)6(.) + a(h^)b{») 

3a{i!))b{f)f)») + 6a{f)f))b{i!)»)b{») + a{f)f)i!))b{»)^ 


3 


Table 3: Substitution law (0 roots) 
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B. Coalgebras, bialgebras and Hopf algebras 


See also |Swe69| | for an introduction to coalgebras, bialgebras, Hopf algebras. 

Throughout this article, k is either M or C, and all tensor products iX> 
are taken over k. An associative algebra over k is a vector space A over 
k equipped with an associative linear multiplication map gi: A ® A —A, 
gL\ a ® b ^ ab. It is unital if, in addition, there is a linear map n: k —>■ A, 
called the unit such that for all /3 G k, a € A, u{/3)a = au{f3) = /3a. 

A is graded if A = 0^o A„ as a vector space and fi{An, Am) C A^+m- 

A coassociative coalgebra over k is a vector space C equipped with a co¬ 
product A:C^C^C. In the notation of Sweedler, we write A(c) = 

E(c) C(1)«)C(2). 

The coproduct is linear and co-associative, that is 
(A (g) /) o A = (/ 0 A) o A, 
where I is the identity map on C, or equivalently 

^(C(l)) ® C(2) = C(1) 0 A(C(2)). 

(c) (c) 

A coalgebra is counital if, in addition, there is a linear map e: (7 ^ k, called 
the counit such that for all c € C, X](c) ^('^(i))c(2) = c = Z](c) '^(i)^(c( 2 ))' 

The (co)algebras considered in this paper are all both (co)associative and 
(co)unital, and we will omit the qualifiers. 

A coalgebra C is graded if C = ©“=0 Cn as a vector space and A{Cn) C 
Cl- A coalgebra is connected if Cq ^ k. 

For V, W algebras, respective coalgebras, F 0 IF is again an algebra, re¬ 
spective coalgebra in a canonical way: 

(ui 0 Wl){v2 0 W2) = V1V2 0 W1W2 


or 

A{v 0 re) = EE U(i) 0 0 U(2) 0 tC(2) 

{v) (w) 

A bialgebra H is simultaneously an algebra and coalgebra such that the 
maps H and u: k —)■ H are morphisms of counital coalgebras 

and A:H H ® H and e: H —)■ k are morphisms of unital algebras, that is, 
for all h,k ^ BH, 


A{hk) = EE ^(1)^(1) ® h(2)k{2) 

{h) (fc) 

A(u(l)) = u(l) 0 u(l) 
e{hk) = e{h)e{k) 
e{u{l)) = 1 
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B is graded \i B = 0^o Bn is graded as an algebra and coalgebra, and 
connected if Bq ~ k. 

For C coalgebra and A algebra, Hom]ii(C', A) equipped with the convolution 
product f * g = g-A ° {f ® g) ° is an algebra. The unit of the convolution 
product is ua^ ec- 

For B a bialgebra, Hom]k(i3, S) equipped with the convolution product is 
an algebra. If the convolution inverse of the identity map exists, it is called 
the antipode, and in this case the bialgebra is called a Hopf algebra. 


Theorem B.l. \Swe6!i. Exercise 1, p. 238] A graded, 


a graded Hopf algebra. 


connected bialgebra is 
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